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Abstract 

We compare the predictions of hybrid inflationary models that produce both adiabatic fluctua¬ 
tions and topological defects to first year WMAP results. We use a Markov Chain Monte Carlo 
method to constrain the contribution of cosmic strings and textures to the CMB angular power 
spectrum. Marginalizing this contribution over the cosmological parameters of a power law flat 
ACDM model, we place a 95% upper limit of 23% on the topological defects contribution to density 
fluctuations, the maximum likelihood being of order 4%. This corresponds to an upper limit on the 
string scale of Gy < 3.2 x 10“^. We also explore the degeneracies between the defects contribution 
and other cosmological parameters. 
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I. INTRODUCTION 


In the 1980’s and 1990’s, cosmologists focused on two plausible scenarios, each of which 
produces scale invariant density fluctuations. In the inflationary model, quantum fluctua¬ 
tions are amplified during inflation and produce adiabatic, gaussian and nearly scale invari¬ 
ant fluctuations [l|. This model is successful in resolving some issues of the Hot Big Bang 
Model and is in good agreement with observations over a wide range of scales . 

The alternative defects scenario rests on the realization that since the Universe has 
steadily cooled down since the Planck time, spontaneous breaking of symmetries must have 
occurred. Each symmetry breaking may lead to the creation of topological defects such as 
monopoles, textures, domain walls, or cosmic strings j3|. A non ne^gible contribution of 
monopoles and domain walls is ruled out by some basic observations [^. However, textures, 
which appear when a non-Abelian symmetry is spontaneously and completely broken, and 
cosmic strings, due to a U{1) symmetry breaking phase transition, are both liable. Even if 
they are not the dominant source of fluctuations, they may still make detectable contribu¬ 
tions to the CMB and large scale structures formation. 

Hybrid inflationary models, e.g., D- and F-term inflationary models, predict that there 
should be both quantum fluctuations, occurring during the inflationary phase, and topolog¬ 
ical defects, created at the end of inflation [^, which also induce anisotropies in the CMB 
by gravitational effects J3. 

Some recent papers [8| on these models conclude that cosmic strings should be responsible 
for at least 50%, and at most 85%, of the amplitude of fluctuations in the CMB anisotropies 
power spectrum. However, other works j9| give a much smaller contribution of cosmic 
strings for the same models, and the contribution of topological defects is highly constrained 
US (best £t with a 18% contribution) from measurements previous to WMAP, in favor of 
adiabatic fluctuations. 

In this Letter, we study the constraints coming from hrst year WMAP results. WMAP 
has indeed given unprecedented precise results on fluctuations in the CMB. It is therefore 
worth studying how models involving topological defects can fit its data, and what are 
the foldings of the constraints hereby established on existing inflationary models, previous 
constraints, and cosmological parameters. 

II. MODEL AND METHOD 


As a non negligible contribution of monopoles and domain walls is ruled out, we only 
consider cosmic strings (CS) and textures (TX) in this work and we assume a model, in 
which the power spectrum of the CMB anisotropies is written; 




model 


= (i-/j-7)c/° + /3c“- + 7q 


■yCS 


-iTX 


( 1 ) 


where (/3,7) G [0; 1]^, /? + 7 < 1 and AD stands for adiabatic. The theoretical power 
spectra of temperature, and electric type^larization for strings and textures have already 


been computed by Seljak, Pen & Turok [ll|. Many tests have been performed on the 
procedure leading to these results, and other independent works Q have given similar 


conclusions, all of them being consistent to within about 10%. These results suggest that, 
for i < 1000, the amplitude of fluctuations in the power spectra of electric type polarization 
and magnetic type polarization for cosmic strings and textures, can be safely neglected 
compared to the one in the power spectrum of temperature. Moreover, the predicted power 
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FIG. 1; CMB angular power spectrum induced by cosmic strings and textures alone (green) and by 
a simple power law flat ACDM model corresponding to WMAP best fit (red). We normalize these 
spectra so that the matter power spectrum normalization for our model is the same as WMAP 
best fit for a flat ACDM model without topological defects. It should be noticed that the points 
of the greeii plot have not been recomputed but directly extracted from the study of Seljak, Pen 
and Turok [ll| . 


spectra for cosmic strings and textures are so close for these values of £ that it is reasonable 
to consider that they are the same. Thus, we assume that and Cj^ are the one shown 
on Fig. m and we can simplify equation ([T]) to 

(J^odei = (1 _ a) + a G™ (2) 

where Cj^ can be as well as Cj^, and 0 < a < 1 is the contribution of topological 
defects to the CMB anisotropies power spectrum. 

To generate the adiabatic spectrum, we use CMBwarp a fast CMB code that agrees 
well with the more accurate CMBFast results. In the parameter range of interest, the largest 
differences between these two codes are generally less than 1%. For the calculations in this 
paper, we assume a simple power law flat ACDM model. We choose to normalize both 
and Cj^ so that, whatever the values of a, the matter power spectrum normalization is 
the one corresponding to WMAP best fit as computed by CMBwarp with WMAP best fit 
cosmological parameters j^. 

Thus, our model is chosen as a flat ACDM model to which we add the contribution of 
topological defects, so that it can be described with only 7 parameters: the baryon density 
(cuft), the matter density (wm), Hubble constant (h), the perturbation normalization (A), the 
optical depth (r), the spectral index (n*), and the contribution of topological defects to the 
power spectrum (a). The idea is to estimate the values of each of these seven parameters 
with 1- 2- and 3-cr confidence contours in the 7-D parameter space. 

To investigate the likelihood space, we use a Markov Chain Monte Carlo (MCMC) method 
[T^ . which enables us to evaluate the likelihood of ~ 5 x 10® models in approximatively 7 
hours. This method is based on Bayes’ theorem, which states that the probability to have 
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the set of parameters s given the observed Ce is 


Vis\Ci) 


V{Q\s)V{s) 
JV{Ce\s)V{s) ds 


(3) 


where V{C(\s) is the likelihood of observing Ci given the set of parameters s, and 'P{s) 
the prior probability density. This theorem tells us we can figure out the value of 'P{s\C^), 
what we would like to do, from the likelihood of obtaining the observed Ci from our set of 
parameters s, that we can compute with the likelihood code provided by the NASA/WMAP 
Science Team on LAMBDA, provided we specify priors. Postulating the equidistribution of 
ignorance, we assume that each cosmological parameter has the same probability to have 
each value between the following lower and upper bounds: 

0<Ub<l 0<Um<l 
0.5 <h< 1.5 0.5 < A < 2.5 
0 < r < 0.3 0 < n, < 2 


It should be noticed, as we will see later, that these priors have no effect on our results, as 
the MCMC never encounters the priors’ boundaries, at the exception of r. For this latter, 
the prior enables the MCMC not to enter unphysical regions of the parameter space. 

Given these priors and Bayes’ theorem, we can then explore the likelihood surface by 
computing the likelihood of the model obtained after each step, until we find the most 
likely region for our cosmological parameters to be. We use the algorithm of ua for this 
exploration, and the chain is stopped when we can make sure it has properly converged (we 
will study this criterion in details as a check of our results), which means we have reached 
and stayed long enough in the region of the parameters space of highest likelihood, and that 
the parameter space has been widely explored. 

In this work, we consider the situations in which the TT and TE components of the power 
spectrum are affected by topological defects. Considering the TT data only, or the TT and 
TE components of the power spectrum, does not change the results, among which is Fig. |21 


III. RESULTS AND DISCUSSION 


These results show that a is more likely to be found between 0.02 and 0.15 at 1-cr with a 
maximum likelihood around 0.04, but scenarios involving a zero contribution of topological 
defects are authorized at 2-a. It is to be noticed that our priors let a free to take any 
values. Nevertheless, we see that a contribution of more than 29% is completely excluded 
at 3-cr. Moreover, we find three main degeneracies between our parameters: when the value 
of a increases, so do the values of Ub, Ug and h. Each of these parameters increases linearly 
with a. 

This can be easily understand by looking at the respective shapes of Cj^ and jsf. 
Indeed, the spectrum induced by topological defects gives a very low contribution to the 
second peak. Therefore, if we keep the cosmological parameters constant and a increases, 
which means that the adiabatic contribution goes down, the height of the second peak will 
diminished, so that it will become incompatible with WMAP results. As our program looks 
for the most likely set of cosmological parameters compatible with WMAP, it has to modify 
the values of the cosmological parameters, in a way that this discrepancy is canceled. In 
other words, the cosmological parameters should be modified in a way that the height of 
the second peak increases whereas the one of the first peak does not change significantly. 
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FIG. 2: We plot the 1-a (black), 2-a (green) and 2>-a (red) confidence contours for first year WMAP 
data for three parameters: h, ojb and Ug- This indicates that scenarios involving a contribution of 
topological defects to the amplitude of fluctuations in the CMB greater than 29% are rejected at 3- 
cr, whereas a zero-contribution is possible at 2-cj. The 1-a contour implies a contribution between 
2% and 15%, with a maximum likelihood of 4%. Each parameter on this graph is degenerated 
with a. 

That is done by decreasing Qj, (which makes the first peak decreases and the second peak 
increases) and increasing h (which makes the first and second peak increases). We can check 
this behavior of h directly on Fig. |21 To verify that the evolution of Qb is the one predicted, 
we have to consider the graph giving u>b, where we can see that tUb increases with a, but less 
than h? does, which confirm the decrease of Vtb- 

Another degeneracy is the one between a and if a increases, so does n*. Indeed when 
a increases, we replace a part of the contribution of the adiabatic spectrum by a contribution 
of topological defects. But, at large scales, P{k) due to topological defects is more concave 
than P{k) predicted by the adiabatic model, which can be taken as P{k) cx: with 

Ug = 0.99 ± 0.04 [^. Thus, if we increase the contribution of topological defects, we must 
increase the convexity of the adiabatic spectrum to keep a shape compatible with WMAP 
results, that is to say that Ug has to increase, which it does. 

If we use the set of parameters corresponding to the maximum likelihood to compute the 
power spectrum, we find a spectrum in excellent agreement with WMAP best ht: the hrst 
two peaks coincide perfectly. Nevertheless, a small discrepancy can be seen for the third 
peak. This is consistent with the model of the power spectrum induced by topological defects 
which gives a quasi-zero contribution after the second peak. Therefore, if the contribution 
of the adiabatic spectrum decreases, it is not compensated by the added contribution of 
topological defects, which makes the third peak decrease. Despite this feature, the ht is still 
compatible with WMAP results. It should also be noticed that this is true for each value of 
a included in the 1-a and 2-a contours; small discrepancies begin to appear in the S-a area. 

During this work, we became aware of the paper 0 , which uses a Bayesian analysis in 
a three dimensional parameter space to constrain the contribution of cosmic strings with 
hrst year WMAP data. In their work, (Um, h and r are chosen to be WMAP best 
ht values, and they study the likelihood of their model involving cosmic strings in a 3-D 
parameters space, the only non-topological parameter involved in this space being rig. The 
results they And with this method are similar to ours concerning the best ht, but quite 
different as far as contours are concerned, and this can be easily explained. Indeed, in 0 
almost ah cosmological parameters are kept constant, at the exception of Ug, with their 
values corresponding to WMAP best ht, which the authors claim is allowed by the fact that 
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the contribution of topological defects should be relatively small compared to the one of 
quantum fluctuations. This is absolutely true for the best £t which gives a contribution 
of cosmological defects of order 4%, the other cosmological parameters being very close 
to those estimated by the NASA/WMAP Science Team, as we have shown in this Letter. 
Nevertheless, this statement is wrong generally speaking. Indeed, we have seen that letting 
all parameters free to vary in the parameter space, shows strong degeneracies between some 
cosmological parameters (namely /i, ujt,, Ug and a) so that it is not possible to choose a 
parameter independently from the others. Doing so necessarily leads to wrong estimations 
of contours, which explain the differences between our results and those given in 0. 

As we were ready to submit this Letter, 0 , in which a work similar to ours is conducted, 
was posted and gives results consistent with our work. However, in this paper, the authors 
use the power spectrum induced by local cosmic strings, instead of the global strings we 
consider. Its shape corresponds to the one we display on Fig. [TJ but the peak that is located 
around i ~ 150 on our power spectrum is to be found at ^ ~ 450 in [l8|. The interesting 
point is that, even with this difference, we still get similar results. This suggests that the 
results obtained by the method we both use is qualitatively insensitive to the model of 
cosmic strings we consider. 


IV. CHECKS 


There are two main risks when using a MCMC. First of all, we have to make sure that 
the chain converges properly, that is to say that we have reached a situation in which the 
likelihood oscillates with a small amplitude around a mean value which will be the center 
of the most likely region for the cosmological parameters to be. Moreover, as the chain is 
necessarily finished, we cannot have access to all areas of the 7-D parameter space. But we 
have to be sure that we sufficiently explore this space, in other words, that we do not forget 
the exploration of entire areas of the parameter space. 

To check these two aspects, we use the approach described in 0 and use the method of 
fl^ . which states that the MCMC gives reliable results if the quantity 


R = 




(4) 


where W is the within-chain variance and B the between-chain variance, is smaller than 1.2 
when computed on ric chains each of which contains 2ne elements, from which we consider 
only the last Ug. In this work, we are more conservative and we have chosen to require 
R ~ 1.1. This result is obtained with 10 chains of 100,000 steps. We also discard the 
first 5,000 points of each chain, to eliminate the so-called burn-in zone, during which the 
stationary distribution might not be reached. The results are absolutely not sensitive to 
this last choice, which is an additional clue showing that we consider enough points in the 
stationary distribution. 


V. CONCLUSION 

The hrst year data of WMAP enable to reject at 3-a scenarios involving a contribution 
of more than 29% of topological defects, and so, in particular, of cosmic strings. In other 
words, D- and F-term inflationary models involving a contribution of more than 50% do 
not match WMAP data, and so, cannot be acceptable scenarios. But most D- and F-term 
inflationary models predict more reasonable contributions of defects. For example, Urrestilla, 
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Achucarro & Davis 0 present a superstring-inspired D-terni inflation scenario which does 
not lead to the formation of cosmic strings. Another example is Rocher & Sakellariadou 

who show that in the context of global supersymmetry, for F-term inflation, and local 
supersymmetry, for D-term inflation, these models can lead to contributions compatible with 
our study. 

We have also shown in this Letter that we cannot reject models involving a low con¬ 
tribution (between 2% and 15%) of topological defects, the best fit of WMAP data being 
given by a 4% contribution of topological defects, which gives a result as acceptable as a 
simple ACDM model. Concerning this point, this work goes in the same way as the results 
obtained by Bouchet et al. US with BOOMERanG data, despite the fact that our best 
fit gives a much smaller contribution for cosmic strings thanks to the precision of WMAP 
measurements. 

Thus, looking at large scales is not sufficient to rule out or confirm models involving 
topological defects for the moment. But as WMAP will soon get more precise measurements 
of the CMB anisotropies for the third peak, it will be possible to increase the constraints 
on a. It should also be worth looking at what topological defects imply at small scales. 
Experimental data usable for these scales will indeed soon be provided by some space or 
ground based experiments such as Planck or ACT, so that we will be able to compare 
theoretical predictions with observations. 
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